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Abstract 

We study numerically stick slip motions in a model of blocks and 
springs being pulled slowly. The sliding friction is assumed to change 
dynamically with a state variable. The transition from steady sliding to 
stick-slip is subcritical in a single block and spring system. However, we 
find that the transition is continuous in a long chain of blocks and springs. 
The size distribution of stick-slip motions exhibits a power law at the crit- 
ical point. 
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1 Introduction 

Stick-slip behaviors in frictional motion have been observed in wide variety of 
systems such as atomically thin lubricant films Q, granular materials || and 
earthquake faults § . Batista and Carlson proposed a phenomenological model 
of one block and one spring, in which the transition from steady sliding to 
stick slip occurs via a subcritical Hopf bifurcation. The block is subject to a 
friction force, which depends on the state of the sliding surface and the velocity 
of the block. A more complicated equation was proposed for the transition 
by Hayakawa. H A chain of blocks and springs was studied numerically and 
experimentally by Burridge and Knopoff for the earthquake faults. || Each block 
is connected to its neighbors with a harmonic spring and the blocks are further 
pulled individually forward through elastic coupling. Each block is subject to 
friction force, which depends only on the velocity of the block. Carlson and 
Langer studied numerically a homogeneous chain of blocks and springs being 
pulled slowly, more in detail. In their model, the sliding friction decreases with 
velocity (the velocity weakening friction) and the steady sliding motion is always 
unstable for all Fourier modes. Jl], || All sizes of stick-slip events are observed 
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and the size distribution exhibits a power law in small scales, but, great events 
occur more frequently than is expected from the power law. 

In this paper, we study a modified version of the single block and spring 
system of Batista and Carlson more in detail. The friction force depends on 
the velocity and state variable, and the friction force increases with the velocity. 
Then, we study numerically a chain of blocks and springs, each of which obeys 
the modified model of Batista and Carlson. In contrast to the Carlson-Langer 
model, in which all blocks are pulled uniformly, only the first block is forced to 
move with velocity V and the other blocks are coupled only with the neighbors 
by harmonic springs. This type of forcing was treated also in the original paper 
by Burridge and Knopoff . || . Elmer studied this type of train model with ve- 
locity strengthening friction and found a size distribution with complex scaling 
exponents.^] Using this block-spring system, we find that a dynamical phase 
transition occurs from the steady sliding to stick slips, and the size distribution 
of the stick-slip events exhibits a power law in all scales at the critical point, 
however, great events occur more frequently as the parameter is decreased from 
the critical point. 

2 Time Evolution and the Period of a Single 
Block and Spring System 

The model equation of a single block and spring system is written as 

K{Vt- x) ~ F(9,v), 

9(l-9)-a9v, (1) 

where M, x is the mass and the position of the block, v — dx/dt is the velocity, 
K is the spring constant, and V is the pulling velocity. The state of the sliding 
surface is expressed by 9. The value 9 = 1 represents that the surface is in a 
solid state, and lower values of 9 imply that the surface is partially melted. 
The term a9v > implies that the value of 9 is decreased, i.e., the surface state 
becomes more melted, as the block velocity v is increased. The sliding friction 
force F(9, v) is assumed to be 

F(6,v) = a6 + j3v. 

The sliding friction contains a standard term (3v which increases in proportion 
to the velocity, and the term a9 which depends on the state variable. There 
are many parameters, however, eq. (1) can be rewritten by scale changes: x — 
a/Mx, v = a/Mi, V = a/MV, K = K/M, (3 = [3/M and a = aa/M as 

^ = K(Vt -i)-6-0v, 

— =9(l-9)-a9v. 
dt 
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Figure 1: Time evolution of z in a single block and spring system eq. (1) at 
M = 1, a = 1, j3 = 1, a = 1, K = 10, V = 0.001 and F s = 3. 



That is, the parameters M and a can be fixed to be M = a = 1 in eq. (1), 
and only the parameters a and (3 are important parameters. We have further 
fixed as a — 1 and studied the model system by changing (3 in this paper for 
the sake of simplicity. We have tried several other parameter values of a, and 
observed similar results. However, the detailed parameter survey is left to the 
future study. When the block velocity is decreased to zero, the block stops. The 
block does not move until the spring force K(Vt — x) goes over a limit value F s 
of the static friction. We have assumed that F s is constant and it is definitely 
larger than o9 in contrast to the model of Batista and Carlson. 

A steady sliding solution is expressed as v = vo = V, 6 = do = 1 — V and 
x = xq = Vt — {1 — (1 — /3)V}/K. The stability can be investigated from the 
linearized equation of eq. (1) 

-KSx -60- (3Sv, 

(1 - 26» - v )Se - 6 a Sv. (2) 

For example, the Hopf bifurcation occurs at V c = 0.106 for K — 10 and (3 = 0.08. 
The steady sliding solution becomes unstable and only the stick-slip motion 
appears for V < V c . However, the bifurcation is subcritical and the stick-slip 
motion appears in wide range of parameters. For example, the steady sliding 
solution is linearly stable for any V > for K = 10 and (3 = 1. Stick-slip 
motion appears at V = 0.001 as shown in Fig. 1. The numerical simulation was 
performed with the Heun method with time step 0.005. The initial condition 
is x(0) = 0, v(Q) = and 8(0) = 0.1. The inverse transition from stick slip to 
steady sliding occurs at V = 0.345 for (3 = 1 and F s = 3, as V is slowly increased. 
The time evolution of v = dx/dt and 9(t) in the slip process hardly depends on 
the pulling velocity V, when V is sufficiently small, which was pointed out in a 
granular experiment by Nasuno et al. H . The time evolution in the slip process 
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Figure 2: (a) Time evolution of v = dx/dt in a slip process for a single block 
and spring system at M = 1, a = 1, j3 = 1, a = 1, K = 10, V = 0.001 and 
F s = 3. (solid curve) and the time evolution by eq. (3) (dashed curve). The 
two curves are overlapped and the difference is invisible in this plot, (b) Phase 
portraits of (v(t),6(t)) by eq. (1) (solid curve) and by eq. (3) (dashed curve). 
Two curves are also well overlapped. 



in the limit of V = is expected to obey 
d 2 x' 

^ = F s -Kx'-6-Pv, 
d6 

— = 9(1-6)- 6v, for t c <t< t e (3) 

where t c and t e are the starting and ending time of a slip process, the position 
variable is changed as x'(t) = x — x(t c ), and the pulling velocity V is neglected, 
and the fact that the spring force K(Vt — x) is equal to F s at the starting time 
t = t c is used. The initial condition for eq. (3) is x'(t c ) — 0, v(t c ) = dx'/dt = 
dx/dt = and 6(t c ) = 1, since the state variable is recovered to 1 during a long 
stick interval. Figure 2(a) compares two time evolutions of v(t — t c ) for eq. (1) at 
V = 0.001 and eq. (3) for F s = 3, and Fig. 2(b) compares two phase portraits of 
{v(t-t c ),6(t-t c )) for eq. (1) at V = 0.001 and eq. (3) for F s = 3. The difference 
cannot be seen in these plots. The interval of the stick depends strongly on V. 
According to the numerical simulation of eq. (3), the time when the block stop 
again is t e = t c + 1.035, and F s — Kx' = —0.432 at the time t = t e . The spring 
force / is slowly increased as / = K(Vt — x) in the stick interval, until the limit 
/ = F s = 3 of the static friction. The interval of the stick process is therefore 
evaluated as At = 3A32/(KV). The total period T is therefore evaluated as 
T = t e — t c + At = 1.035 + 3.432/ (KV). Figure 3 displays numerically obtained 
values of the periods of the stick slip motions and the theoretical curve. These 
numerical results support the model equation eq. (3) for V = 0. We have shown 
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Figure 3: Period of the stick slip oscillation as a function of the pulling velocity 
V. The marks denote numerical results and the dashed curve is T = 1.035 + 
3.432/ (K V). 



some numerical results for specific parameter values, however, the qualitative 
properties do not depend on the parameter values. 

3 Dynamic Phase Transition in a Long Chain of 
Blocks and Springs 

Next we consider a chain of N blocks and N springs. The model equation is 
written as 

K(xi+i - 2xi + Xi-i) - F(8i,Vi), 

8,(1 - Oi) - + D(0 i+1 - 28, + 9i-!), i = 1, 2, • • • N (4) 

where Xi, Vi and 6% represent respectively position, velocity and state variable 
for the ith block, K is the spring constant, and a diffusion type coupling with 
strength D is assumed in the equation of the state variable. The sliding friction 
is assumed to be F(8 i} Vi) = 0, + 0Vi for nonzero v,. Once the ith block stops, 
the block does not move until the spring force goes over the maximum value 
F s of the static friction as the case of the single block-spring system. The 
first block is assumed to be pulled with a constant velocity as x = Vt. The 
steady sliding solution is written as Vi = vq = V, 8, = 6q = 1 — V and Xi = 
Vt -'{1 - (1 - 0)V}/K ■ {N 2 - (i - N) 2 }/2. The steadily sliding state is 
nonuniformly strained. The stability of the steadily sliding solution can be 
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Figure 4: Real parts of the largest two eigenvalues for the Fourier modes with 
wavenumber k for eq. (4) at (3 = 1.5, 1 and 0.5. 

investigated from the linearized equation of cq. (4) 

K(5x i+ i - 2Sxi + Sxi-i) - 58i - f35vi, 

(1 - 29 a - vo)59i - e Svi + D(S6 i+1 - 259, + 59^). (5) 

The real parts of the largest two eigenvalues for the Fourier modes with wavenum- 
ber k are shown at three parameter values 0.5,1 and 1.5 of (3 for K = 10 
and D = 5 in Fig. 4. The steady sliding solution becomes unstable against 
long wavelength perturbations for (3 < 1. The parameter (3 = 1 is a crit- 
ical value. We have performed numerical simulations with TV = 399 for 
K = 10, D = 5, V = 0.001 and F s = 3. No-flux boundary conditions are 
assumed: xjv + i = xn, #at+i = 9n and 9\ = 02- The initial conditions are 
Xi(0) — 0, Vi(0) = 0.1 and 6^(0) = 0.1. Irregular stick-slip motions are observed 
initially. A steady sliding state is finally obtained for (3 > 1 after a long run, 
however, stick-slip motions are maintained for (3 < 1. Figure 5(a) compares a 
snapshot pattern of Xi{t) at t = 2 x 10 7 for j3 — 1 with the steady sliding solution 
Xi = Vt - {1 - (1 - /3)V}/K ■ {N 2 -{i- N) 2 }/2. The good agreement implies 
that the stick-slip state is very close to the steady sliding solution at /3 = 1. 
The average value of the total length L = {x\ — xm) is displayed in Fig. 5(b). 
For (3 > 1, the steady sliding solution appears and the total length is close to 
L = {1 — (1 — (3)V}N 2 /2. In the stick slip phase, the total length decreases 
continuously from the value of the steady sliding solution. These results sug- 
gest that the transition from the steady sliding phase to the stick slip phase is 
a continuous transition. Chaotic motion appears at the onset of the instability. 
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Figure 5: (a) Snapshot pattern of Xj at t = 2 ■ 10 7 (solid curve) and the steady 
sliding solution (dashed curve) for (3=1. (b) Average total length L = {x\— xn) 
as a function of (3. The dashed line denotes the total length for the steady sliding 
solution. 



The slip motion starts at i = 1 in most cases and the rupture process propa- 
gates and ends at a certain position. There are a long interval of stick between 
two slip motions for small pulling velocity V. The slip events occur in a variety 
of sizes. In small slips, the rupture process ends at small i. In large slips, the 
rupture reaches i = N. Figure 6 displays time evolutions of a large-size slip near 
t = 19100 and an intermediate-size slip near t = 40400, where the displacement 
Axi = Xi(t + At) — Xi(t) during a time interval At = 20 is plotted. The leading 
edge of the rupture propagates with a nearly constant velocity c = 1.43 for 
large i, which is fairly smaller than the sound velocity y/~K = 3.16. (Langer and 
Tang studied rupture propagation in their homogeneous model, epicenters are 
randomly distributed, and the rupture propagates in both directions from the 
epicenter. Jl0| ) 

Figure 7 (a) and (b) display a distribution of the event size and the time 
evolution of the event at /3 = 1. Here the event size is defined as the total 
displacement S = Yli=i during each slip event. For the critical parameter 
(3=1, the distribution exhibits a power law with exponent nearly 1.5. Vari- 
ous size of events occur as is seen in Fig. 7(b). The largest size seems to be 
determined by the system size. Figure 7 (a) and (b) display a distribution of 
the event size and the time evolution of the event at (3 = 0.2. For (3 = 0.2, 
the distribution seems to obey a power law with exponent about 1.6 for small 
events, however, the distribution deviates from the power law and very large 
events occur more frequently than is expected from the power law. Very large 
events with similar sizes occur frequently. 
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Figure 6: (a)Time evolution of the displacement Axi = Xi(t + At) — Xi(t) in 
a time interval At = 20 in a large slip event for — 1. (b)Time evolution of 
the displacement Ai, = xi(t + At) — Xi(t) in a time interval At = 20 in an 
intermediate slip event for f3 = 1. 



4 Summary and Discussion 

We have proposed a model which exhibits a phase transition from steady sliding 
to stick slip. For a single block spring system, we have reconfirmed that the 
time evolution in the slip process hardly depends on the pulling velocity V when 
V is sufficiently small. The equation of the time evolution and the period have 
been evaluated in the limit V = 0. Hayakawa evaluated similar quantities for 
his model through the reduction to the Amontons-Cloulomb model in the limit 
of V = 0, neglecting the friction force ||. We have evaluated the period in the 
limit V = 0, retaining the friction force. In a long chain of blocks and springs, 
the blocks are connected to the neighbors only with harmonic springs, and only 
the first block is pulled with a constant velocity V. We have found a continuous 
transition from steady sliding to stick slip. In most cases, slip processes, i.e., 
ruptures start from the first block. At the critical point, the distribution of slip 
size obeys a power law, however, large events with a typical size occur more 
frequently than is expected from the power law for small (3, which is consistent 
with the model of Carlson and Langer. Our model is a deterministic system, 
however, it exhibits a dynamical phase transition similar to the thermodynamic 
phase transition including the critical phenomena. Bifurcations and chaotic 
behaviors were studied in a model of two blocks and the connecting springs 
by Vieira pi] , but our model of a long chain of blocks and springs exhibits a 
continuous transition from a stationary sliding state to a spatio-temporal chaotic 
state. Our model suggests that the size distribution depends on a parameter, 
and the power law distribution is observed near the transition point, although 
direct relevance to physical systems such as earthquakes is not clear yet. 
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Figure 7: (a) Histogram of size distribution of the total displacement S = 
YliLi A-Xi during each slip event at /3 = 1 in logarithmic scales, (b) The total 
displacement S as a function of the event time at (3 = 1 
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Figure 8: (a) Histogram of size distribution of the total displacement S = 
YliLi during each slip event at (3 = 0.2 in logarithmic scales, (b) The total 
displacement S as a function of the event time at j3 = 0.2 
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